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1. Introduction 

Ions in an ion trap usually behave as a gas: they try to fill all the available volume, 
defined by the trap's electromagnetic fields. It was only in storage rings where the 
ions could be forced into well localized bunches. But even there the ions' mutual 
repulsion and the inevitable spread in their kinetic energies tend to distribute the ions 
throughout the whole ring. 

Recently a surprising cooperative behavior of the ions was discovered in a Zajfman 
trap , an ion trap resonator built of only two opposite electrostatic mirrors and two 
focussing lenses [||: the ions are injected into the trap from an ion beam and this 
injected bunch of ions bounces back and forth between the two mirrors but docs not 
diffuse in the trap. The ions, which all have the same charge and also slightly different 
energies and trajectories through the trap, seem to "stick together"; they are not 
only confined to the trap volume, but also seem to be trapped within the bunch, 
which oscillates back and forth like a single macroscopic "super particle" . This bunch 
works as a "dynamic inner trap" , itself moving between the static electric fields of the 
external trap. 

As a first application of this "Zajfman trap" operating in the bunching regime a 
high resolution Fourier transform mass spectrometer has been demonstrated with a 
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resolution of Am/m w 7 x 10~® ||^. This value was until now only achieved in storage 
rings This high resolution is made possible by observing the oscillations of the 
macroscopic bunch through the trap over time intervals of several tenths of a second. 
The observation time is limited by the trapping time of the ions in the inner trap (the 
bunch) , which is determined by the number of collisions between the ions and residual 
gas atoms in the vacuum chamber Q . 

The operating regime of the trap and the stability conditions, under which this 
counter intuitive "self-bunching" , or synchronization, occurs, were first explained by 
Strasser et al by relating the ion dynamics to the so called negative mass instability 
[^: due to the special dispersion of the trap the period of an ion increases with its 
energy. This can be described by a negative effective mass of the ions. Consequently 
the ions have to repel each other in order to synchronize their motion. 

A microscopic explanation of the observed synchronization for two identical ions 
was given by the authors in reference which we refer to as "part I" in the following. 
There the problem is described by a set of simple mappings for the evolution of the 
ions' relative coordinate through the different parts of the trap. The ions' interaction is 
incorporated phenomenologically as a time delay. With this approach we confirm not 
only the already known stability conditions but describe the microscopic process taking 
place: synchronization is a continued alternation of energy and position exchange 
between the two ions. 

In this paper we show how the model of two identical ions can be generalized 
to describe the dynamics of the whole bunch. For this we associate one of the ions 
with a "test ion", which moves through the cloud of the other ions, and the other 
describes the motion of the bunch's CM. With the two identical ions we did not need 
to specify the interaction to explain the mechanism, but now the interaction potential 
that the test ion feels is an important part of the effect: it is derived from the charge 
distribution in the bunch, i.e., the test ion's average position in the bunch, which in 
turn is determined selfconsistently from the motion of the test ion inside this bunch. 
This test ion case is a specialization as well as a generalization of the two identical ions 
description: by specifying the interaction we can describe the generalized behavior of 
arbitrarily many ions through a mean field treatment. 

The paper is organized as follows: Section ^ reviews the mapping approach of 
part I. We show that the approach can be extended from the two-ion case to N ions. 
The key issue left unresolved in Section |2] is how to calculate the time delay of the 
test ion in the bunch. This question is taken up in Section |^, where we show how 
this parameter can be obtained from a mean field treatment of the bunch potential. 
A central result of Section || is the derivation of the spherical geometry of the bunch 
at the turning points. In section ^ the dependence of the bunch length on the various 
parameters of the trap and the ions is derived. In section |^ the momentum spread of 
the bunch is calculated. Then, in section ^ we show how the experimental criteria for 
bunching deduced from the observations by Pedersen et al are connected to our 
explanations. The results are summarized and further developments are sketched in 
section |^. 
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2. Review of the mapping model 

2.1. Coordinate system, transit times and the dispersion relation 

In this section we briefly review the basics of the two ion case. For information about 
the trap and the experimental findings the reader is referred to references |0| and 
For more details about the mapping model we refer the reader to "part I" 
Consider two ions with the simplified trap potential: 

f when \x\ < ^ 

Vix) = \ ^ ] (1) 

\ (|x| - -DF when |x| > f. 

F is the constant gradient of the mirror fields, L is the field free distance between the 
two mirrors. 

In general the two ions have masses mi and m2 and charges gi and q2 , respectively. 
As synchronization shows up in the distance between the ions, we use center of mass 
(CM) and relative coordinates. With the total mass M = nii + m2 and the reduced 
mass u = "''-."'^ we define 

mi+m2 

U p if (2) 

X = Xi- X2 J y X2= H - -jf-X . 

The corresponding momenta are 

P=Pl+P2 I f = 

^'=^P1-^P2 j \ P2^'ifP-p. ^ ^ 

R and P are the CM coordinate and momentum, respectively; the relative motion is 
described by x and p. The ions' positions and momenta arc denoted by xi, pi and X2, 

With the ion-ion interaction W{x) the Hamiltonian for these two ions reads 

2 2 

H=^ + ^+ qiV{xi) + q2V{x2) + qiq2W{xi - X2) 
zmi 2m2 

= ^ + + qiViR + '^x) + q2V{R - ^x) + q^q2W{x) . (4) 

To separate CM and relative coordinates when both ions arc in the same part of the 
trap potential we need that qiV{xi)+q2V{x2) = (qi +q2)V{R). With the potential of 
equation we get the condition that both ions must have the same charge to mass 
ratio: 

gim2 - 92^1 qi q2 

rr ^ = (5) 

M mi m2 

For two identical ions this condition is trivially fulfilled. Here we describe a bunch of 
A'' + 1 identical ions; then the above condition is also fulfilled, if we take one of the 
ions to be a test ion of mass m and charge q and identify the other with the remaining 
> 1 identical ions of the bunch: 

mi ~ m m2 = A^m (6) 
qi^q q2 = Nq (7) 
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The number of ions N is consequently a positive integer, which can be iV = 1, too. 
Now, with M = (N+l)m and /i = the CM and relative coordinates (equations 



and M)) become: 



A^+l -^1-^+^+1^ 
X = Xl — X2 J 1 .T2 = i? — jv+X^ 



and 



P = Pi+P2 \\pi=nTiP+P 



^ ^ T • (9) 



In part I we combined the geometry of the trap and the ions' mass and energy 
into two parameters, which are derived from the times that the ions spend in the 
mirror, i.e., at |x| > L/2, and in the central part of the trap, respectively. 

An ion with the laboratory energy Eq and the momentum pq = \/2mEo spends 
the time Tm inside the mirror: 

^" ^ qF~ NqF ^^^> 
Tm is one of the two parameters mentioned above. 

When the ion turns around at T',„/2 it has penetrated into the mirror potential 
for the distance 

X,n = (11) 

qF 

The ions have the velocity po/m « P/M and therefore need the time 

Lm LM , , 

Tf = = := aTm 12) 

to pass through the central field free region of the trap. 

It turns out that not the times themselves but their ratio a — the second 
parameter — determines, if synchronization occurs. When we express a in terms 
of the experimental parameters, where U is the acceleration voltage, i.e. Eq = qU, we 
get: 

LqF LF 

a = = , (13) 

AEo 4[/' ^ ' 

With T„j and Tf we calculate the total time T for one period as 

T = 2Tm + 2Tf = (14) 

(15) 

This form highlights the importance of the bunching effect for mass spectrometry 
applications: the period of the bunch, i.e., the synchronized ions, which can be 
measured very precisely for a huge number of revolutions, is directly related to the 
ratio ^. 

The factor a also shows up in the dispersion of the trap, i.e., in the dependence 
of the period of one ion on its momentum. With equations (jl^) and ( |l3| ) the dispersion 
is a simple relation between the "geometric" parameters: 

|..l(l-„) = i(«„.-L) (16) 

For (1 — a) > the period T increases with the ion's momentum po = \/2mE() and 
vice versa. 
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2.2. Independent ions: The momentum kick and the basic mapping structure 

In our model potential the CM and the relative motion separate for nearly all the time, 
except for a short interval . This is the time between when the first ion crosses the 
kink between the flat part and one of the mirrors at ±-|- and when the second one 
passes that point, too. During that time the ions arc accelerated relative to each other 
by the trap potential. For the case of ions going into the mirror at L/2 it is: 

rk^^ (17) 

After this time the momentum p2 of the second ion will be the same, P2 = P2j but pi 
is slowed down to 

p'i=Pi+Tfe =pi qF. 18 

V OXi ) P2 

The two ions have almost the same velocity, so — is approximated by — . The relative 
momentum p = j^jjPi — j^^P2 now changes with equation ( pl)| ) as 

NmqF 2/i 
V ^P- . -I — x=p-—x. 19 
N + 1 Po T„i 

As Tfc is short, the relative distance x does not change during this time in our 
approximation; also the ions' very weak interaction can be neglected [Q. This sudden 
change of the momentum (|l^) can be formulated as a mapping K, of the relative 
distance and momentum x and p just before this kick onto the values x' and p' right 
after the kick: 

Without the ions' interaction the relative coordinate evolves freely in the central 
part and inside the mirrors, leading to the two mappings J- and M, respectively, which 
differ only in the length of their time interval: 

: )-( : ) = ("^/''') (21) 



The complete dynamics of two noninteracting ions in their CM system can be 
described by these three mappings /C, J- and M. They consequently form the 
foundation of our mapping model. With the ions' interaction incorporated into this 
picture we are able to determine the stability conditions for bunching and, as we will 
show in this paper, also describe the extension of the whole bunch with respect to the 
external parameters. 



2.3. Including the ion-ion interaction: the time delay and the stability condition 

As shown in part I 0|, the ions' motion is synchronized by their repulsion when the 
ions interact in the mirror. Without further specifying the form of the interaction 
potential we can model its effect by a time delay r„i that affects the relative motion 
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when the ions' paths cross each other. Wc include this delay into the mirror 
mapping as 

Now we can determine the stability of the motion by calculating the eigenvalues of 
the linearized mapping 

■^1/2 = ^1/2 ®K.®M' ®K.® J"i/2 (24) 

of a half period. They have the form 



Ai/2 = -l + 27±2V7'-7, (25) 
where we applied the abbreviation 

7 = e(l - a) with e = ^ (26) 

Equation ( p5| ) was one of the central equations of part I. Analyzing the two ion case 
in this paper it will take on a new significance in describing the stability of the N-ion 
bunch. 

The three mappings IC, M' and T are independent of the ratio of the two masses 
mi and m2; consequently the behavior of the relative coordinate and the stability 
criterion are the same here in the test ion case (toi <C 7712) as for the two identical 
ions, which was discussed extensively in part I 0. 

The eigenvalues A1/2 (|25|) are complex for < 7 < 1. Then the test ion is stably 
synchronized with the bunch |Q — which then will be stable, as the test ion describes 
the motion of all the other ions in the bunch, too. This requires that the dispersion be 
positive, i.e., a < 1, and the delay in the mirror be positive, r„i > 0. For a negative 
delay, < 0, which would lead to synchronization for a > 1, in general an attractive 
interaction is necessary, which does not occur for equally charged ions. We will in the 
following only consider the case of a < 1 and > 0, i.e. the interaction between 
the test ion and the other ions of the bunch is repulsive and it takes place around the 
bunch's turning point in the mirror. 

To explain the observed synchronization in the two ion model of part I |Q it was 
sufficient that the ions experience a time delay, when their paths cross in the mirror. 
The magnitude of this delay did not matter, only its sign. However, for understanding 
the geometry of the N-ion case the magnitude is important as well. In section || we 
will show that the time delay can be obtained once one knows the bunch geometry, 
while the bunch geometry in turn determines the delay. Thus, these two quantities 
are determined self-consistently. 



2-4-. Incorporating the effect of off -axis motion 

There is yet another difference between our one dimensional model and the real trap 
that wc need to include: in the real trap even two ions with exactly the same energy 
may have different periods, as the length of their orbits varies with the radial distance 
off the trap's central axis. Pedcrscn et al characterized this spread of isoenergetic 
periods by a time difference AT; If this spread is too big, synchronization is 
suppressed in the experiment. 

As our one dimensional model does not allow for off-axis orbits, we include 
this spread as an externally given parameter in the following way, explained for two 
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(a) 



• • ■ 
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Figure 1. The phase space positions of the ions in the injected bunch (a) serve 
as starting points for the ions' oscillations around the center of the bunch in the 
stable regime (b): no matter how small the energy (momentum) spread of the ion 
source is (a), the bunch's internal energy in the trap depends on its length. 



identical ions: In the central part T/ is effectively independent of whether a given 
trajectory lies on or off the trap's central axis. So this difference in the period is 
accumulated in the mirrors. There a delay of one of the two ions with respect to 
the other due to the spatial dispersion leads to a bigger separation between the ions, 
which is equal to extending the time that the CM spends inside the mirror by ^5 if 
the total spread after one period is AT^ — 2Ti we consequently have to replace the 
time interval Tm in A4' (^3|) by Tm + This is valid for a test ion and the bunch, 
too. The mirror's mapping then takes on the following form: 

x\ f x' \ / .T + i(T,„-(r„-^))p 



This form reduces to equation ( p3|) when we define an effective delay = t,„ — 
Consequently the stability criterion > for a < 1 has to be replaced by 

AT 

t'^>G ^ ^™>^- (28) 

The stability of the bunch is independent of the time spread ATi, (see caused by 
the finite energy resolution of the ion source: when the two ions are in the synchronized 
mode they oscillate around each other, exchanging momentum all the time. Therefore 
even a beam which had been injected without any energy spread is "heated up" 
corresponding to its length, as explained in reference ||^ and figure g. The energy 
spread of the initial injected bunch can be seen as a snapshot of the different ions 
somewhere on their orbit around the bunch's CM. 

This heating is related to the time delay Tm, which we will determine in the 
following section. When we then have actual numbers for we will calculate the 
momentum and energy spread of the bunch. 

3. Calculating the time delay and the bunch length 



From the explanations of the previous section 2.4 on the off-axis motion we conclude 
that the bunching phenomenon is more stable against perturbations when the time 
delay is larger: the observed bunch consequently will be the configuration that 
maximizes the delay for the given experimental parameters. 

Thus, the key to understanding the properties of the N-ion bunch is to quantify 
the size of the time delay . This time delay depends on the bunch geometry, while 
the bunch geometry and stability in turn depend on the time delay. Thus these two 
quantities have to be determined self-consistently. 
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(i) (ii) (iii) 
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Figure 2. The ions' motion in the mirror is divided into three stages: (i) while the 
bunch penetrates into the mirror the bunch contracts, but the ions do not change 
their relative ordering, (ii) When the bunch has reached its minimal length A at 
the turning point the ions cross over the potential barrier to the other side and 
(iii) on its way out of the mirror the bunch expands again symmetric to (i). The 
behavior of the explicitly denoted test ion is representative for the motion of all 
the other ions. 

To do SO we use the modified two ion model introduced in the previous section: 
one of the ions with a mass mi = to is tlie test ion. It moves relative to the other 
ion with TO2 = Nm, describing the CM of the other N ions of the bunch. The ions' 
interaction, and therewith the time delay, which we did not have to specify in the two 
ion case, is now derived from the bunch geometry. 

3.1. Modified three-stage process 

Consider first the case of non interacting ions. Then we can divide the bunch's behavior 
in the mirror into three distinct steps, see figure ||: (i) In the first stage, the bunch 
enters the mirror. At this stage the ions get a compressing momentum kick, which 
for each ion is proportional to its distance from the CM, see equation Given 
the macroscopic length of the bunch, the momentum kick is much bigger than the 
spread of the ions' momenta in the central part of the trap. Therefore, to a very good 
approximation, the ions do not change their relative order, while the bunch contracts 
to some minimal length, shortly before the CM turns around in the mirror, (ii) In the 
second stage the CM turns around and the ions move through the now (stationary) 
bunch to the other side, (iii) In the third stage the bunch expands again on its way 
out of the mirror, until the second kick at the exit of the mirror undoes the first 
momentum kick and stops the expansion. 

Now we have to incorporate the ions' interaction into this process. When their 
interaction is weak compared to the momentum kick, this three stage process will 
pertain. Since the effect of the ions' interactions scales with their density their 
repulsion is most effective when the density is the highest, i.e., when the bunch 
is compressed to its minimal length at the turning point during the second stage. 
We incorporate the ions' interaction into this stage only and arrive at a modified 
three-stage description: the first and last stage, where the bunch is compressed and 
expanded, respectively, remain unchanged. The modification is to the second stage 
only, during which the ions feel the potential of the compressed bunch and accumulate 
a time delay while they move to the other side of the bunch. Consequently the time 
delay and, therefore, all properties of the stable bunch are determined by its geometry 
at the turning point. 



Synchronization in the "Zajfman trap" II 



9 



3. 2. Transformation between the bunch lengths in the central part of the trap and at 
the turning points 

As discussed in the previous paragraph, the motion of the ions during the second stage 
is determined by the momentum kick that they got at the entrance to the mirror. 
The magnitude of this momentum kick is in turn determined by the ions' respective 
distances from the bunch center when the bunch enters the mirror. We therefore have 
to derive a transformation for the first and third stage relating the bunch geometry in 
the central part of the trap to its extension at the turning point. 

For noninteracting ions the length of the bunch in time is the same both at the 
observed length and at the minimal length A. It is identical to the time interval , 
which defines the momentum kick ([it]), expressed here with the average momentum 

rk = ^ (29) 

The minimal length of the bunch occurs when the CM turns around in the mirror. 
Then the first ion, being in front by r^., has already returned from the turning point 
by the distance 

while the last ion is still the same distance before X,„ (pi]). Here we assume that the 
mirror potential at the turning point is linear over the (small) length 2A of the bunch. 

In the two ion case Q there was no need for some minimal bunch length, as the 
crossing time of two ions is always well defined, but for many ions the spatial extension 
of the continuous bunch translates into a length in time, during which the ions cross 
the bunch. 

With the above form of (KS) and -Bo — ^ we get 

A=^(-^] Ll (31) 

Note that the form of the trap potential between the kink and the turning point docs 
not influence the minimal bunch length: seen in the laboratory frame all ions "climb 
up" the same potential ridge with essentially the same initial energy Eq; therefore 
in our approximation their spacing in time remains constant all the way up to the 
turning point independent of the actual form of the mirror potential. We only require 
a linear potential right around the turning point. 

For the following we nevertheless use the special form of our model potential (|l]) 
and come back to the general form only after we have analyzed the behavior of the 
observed bunch in our model potential. Then equation ( pT| ) becomes: 

A = ^^Ll (32) 

2por„i 

With the values of the initial experiment i.e., Ar+ ions with Eq = 4.2 kcV and 
F = 80 kV/m, a bunch of, e.g., Li, = 2 cm is compressed to A = 0.95 mm. 

With the linear slopes at the turning point the bunch is contracted linearly during 
the first stage. Therefore an arbitrary position x inside the uncompressed bunch of 
length Lf, is transformed into Xi in the contracted bunch, i.e., when it has the length 
A, as: 

A mLh 
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3.3. Calculating the time delay 

According to the three stage description explained above we have to deal with the 
following process: (i) the linear contraction translates an ion from a position x in 
the expanded bunch to the position Xi in the contracted bunch (|33|). (ii) From this 
position the ion starts to cross over the now stationary bunch potential to the other 
side until it reaches —Xi. (iii) Finally, the linear interaction free expansion picks it up 
again. 

Consequently the time delay due to the repulsive bunch potential is defined 
as the difference between the time T that the test ion needs to cross over the bunch 
of minimal length A, starting from an initial position Xi up to the symmetric distance 
—Xi, and Tq, the time needed for the same distance 2xi without any potential: 



Note that with the interaction the spatially symmetric positions Xi and —Xi are 
not reached at times symmetric to all ions are delayed in the central stage, 

consequently they all reach their symmetric position on the other side of the bunch 
later than without interaction. The third, expanding stage only starts when all ions 
have changed to their respective side of the bunch. Of course, this separation into 
three successive stages is an idealization, which will not be observed strictly in the 
experiment. 

To calculate the times T and Tq we need the momentum of the ion at the end of 
the first stage, i.e., at Xi. In addition, to determine T we need the explicit form of the 
bunch potential. These quantities will be calculated in the next two subsections. 

3.3.1. Calculation of the momentum of the test ion To evaluate equation ( p^ ) we 
start with the momentum p{xi): it is the sum of the ion's initial momentum relative 
to the CM in the central part of the trap plus the momentum kick at the position 
X. For a macroscopic bunch the kick is much bigger than the momentum spread: an 
energy spread of 10 eV in the central part at an ion energy of 4.2 keV corresponds to 
a momentum spread of about 5.7 au for Ar^ ions; during the time Tm this momentum 
difference widens the bunch by about 0.25 mm, much less than the distance 2if, of 
a few cm, which is traversed during the same time due to the momentum kick. We 
consequently neglect the momentum spread and calculate the ions' energy relative to 
the CM from the kick alone. 

The momentum change of one ion relative to the bunch at the position x due to 
the kick was shown to be (equation ([l9| ) with ^ = m) 




(34) 




The corresponding energy relative to the bunch is e = = . With the 

transformation x — ^f^J" Xi ( |3^ ) we calculate the energy at the contracted position 
Xi as 




'6 



With e{xi) we can now determine Tq (cf. (^)): 
2xim Lhm 
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Tq is independent of Xi because due to the linear contraction of the bunch the 
momentum kick is proportional to the starting distance Xi. 

3.3.2. Approximating the bunch potential To calculate the time T that the test ion 
needs to pass over the bunch we need the ions' distribution in the bunch. From this 
we calculate the bunch potential and then the delay, which in turn determines if this 
specific bunch is stable. But if it is not stable for all Xi, the ions will redistribute, 
modifying the bunch potential. In a more exact treatment the bunch potential would 
have to be iterated so that for all ions in all parts of the bunch the stability criterion 
is fulfilled. 

To keep the description simple we will not allow for an arbitrary bunch form 
but explicitly make the (mean field) ansatz that the ions are spread in a Gaussian 
distribution around the CM all of the time, i.e., in the central part of the trap and 
during the whole evolution through the mirror. With this assumption we only need 
one parameter to describe the width of the bunch. The charge density p{x) of the 
contracted bunch then will be, with the total number of ions N in the bunch and its 
length A at the turning point: 

P(.^) = ^e-^V2A= (38) 
V2t:A 



The ions interact via a repulsive Coulomb potential jjy. To allow them to pass by 
each other in our one dimensional model this singular potential is replaced by a so 
called "softcore" Coulomb potential ||lT[| : 

^ TO 

The " softcore parameter" or " impact parameter" d is a measure for how close the ions 
have to come when passing each other; it describes the radius of the ion beam in the 
trap potential. 



With this ansatz and the Gaussian density profile (38) we calculate the bunch 
potential Wb as 



We could now insert this potential into (^4|), but then the delay for a specific ion 
depends on where it is located in the bunch. Therefore we concentrate for the following 
on the ions around the center of the bunch. Introducing V{xi) as the quadratic 
expansion of Wb around the center of the bunch, we have 

Wbix,) « y(xO = WbiO) - (41) 



Clearly, we have 



= muj^ , (42) 

x=0 

a relation which will be needed below. 

With this harmonic approximation we get for the energy Ej = E{xi) of the ion 
above the top of the bunch potential: 

= e{x,) + V{x,) - Wb{0) = (^M _ ^2 (43) 
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From this energy we may calculate the position dependent momentum as 



2E 

p{z) = mwW ^ + 2;2 . (44) 

V mcj^ 

Remember that Xi denotes the starting position of the ion at the end of the first 
compressing stage, just as it starts to cross the bunch during the second stage, while 
the integration variable z denotes the (time dependent) position during the second, 
the crossing stage. 

Using equations ( p^ ) and (^) the time T that the ion needs to cross over the 
approximated bunch potential V{xi) (^l|) evaluates as: 



2 

T = -Arsinh 

iU 



(45) 



Here we reinserted the initial energy of the ion Ei ( p3| ) into the time dependent 
momentum p{z) (^). Again, as for To, T, and therefore t„i, is independent of Xi 
due to the harmonic approximation ([H]). This approximation allows us to focus on 
how the delay depends on the parameters of the trap and the bunch without having 
to deal with a whole range of initial conditions Xi. Remember that we already had 
neglected the ions' momentum against the momentum kick. 

The second derivative in mui^ ( ^2|) and the integration in Wi, (^ ) can be 
interchanged to further evaluate the curvature of the bunch potential's top: 



q^N f dz „_,V2A^ (-^ _ 



= 4^ldye-y'-'^ (46) 

By changing the integration variable to y = z/A the integral depends only on the 
fraction 8 = d/A, i.e., the ratio between the beam diameter 2d and the (total) length 
2A of the bunch. In the following we will abbreviate the integral as 

I{d/A) ^ I{6) = [dyc~y'/' ]:/ . (47) 



From the expressions for t„i according to (|3j), (|37|), ( [ISD and (|4(i| ) we could 
now calculate the time delay numerically and learn about its behavior. But 
these expressions are too complex to easily "see" the dependencies on the various 
parameters. We will therefore derive a simplified form thereof in the following section. 



3.3.3. Linearizing Tm In the previous section we derived the equations necessary to 
calculate the time delay numerically for a given trap configuration. However, these 
equations can be simplified so that the central stability criteria become much more 
obvious. 

We can linearize the above set of equations by expanding both the Arsinh and the 
square root in equation (|45|). This approximation is valid for Ltmuj <C 4po, i-C, when 
the ions' energy from the momentum kick is much higher than necessary to cross the 
bunch (cf. equation ([43[)). Then the resulting delay is small. From equations (p8[), 
( pO| ) and ( ^ ) we see that Wb, and therefore muj'^, scales linearly with TV. When N is 
kept small enough this condition is fulfilled even for a short bunch or a thin beam. 
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We first set (1 — x) w 1 + x/2 and then Arsinh(a;) w a; — x'^/Q to approximate 
equation (ES) as: 



T; 



—Arsinh 

LU 








2po V 


2po ) 



(48) 



12 



32p2 

-0(c^^) (49) 



The first summand in the second line equals Tq (^). With equations (|T(i|), (^), ( ^6| ) 
and ( ^T] ) we finally arrive at an expression for r™, which depends only on parameters 
that describe the bunch at the turning point in the mirror: 



3/2 

- I 



12" y^j 

Jmq N 



mu^ (50) 
I{6) (51) 



12V2^^^3/2 

Finally two degrees of freedom are left, which the bunch can adjust to achieve the 
most stable configuration, i.e., the biggest Tm- the number of ions N and the ratio 
6 ^ d/A between the beam width and the bunch length at the turning point. All 
other quantities are fixed in the experiment. 

From ( |5l| ) it is readily seen that the delay is proportional to the number of ions 
TV in the bunch: 

Tm oc N (52) 

The bunch can consequently stabilize itself by taking up more ions; or, put the other 
way, the bunch has the tendency to keep all ions together. It is stable with respect to 
the number of ions. 



3.4- Spherical bunch geometry 

The dependence of r,„ on the beam diameter d and the bunch length A at the turning 
point in the mirror is less obvious, but if there is a "most efficient" value of A, i.e., one 
for which the delay is maximal for a given d, then for this configuration the derivative 
vanishes. Combining all the prefactors which are not related to or ^ into one 

constant C = ^gV^Tf^/a ' bom (^) and (|4^) : 

It will be shown later that the extremum at ^fx- = is in fact a maximum. Since 
C and A always have a finite value, the maximum of Tm is consequently determined 
by that value of 5, for which the integral vanishes. This happens for 5 ~ 1.01689 . . ., 
which is about unity. Consequently, Tm is maximal for Ak d, i.e., when the bunch at 
the turning point has a spherical form. With A ^ d the integral ( [47| ) in (^Tj) takes on 
a constant value of X(l) = 0.5778 . . .. 

This value of S of nearly unity becomes more clear when we artificially replace for 
a moment the softcore interaction in the convolution of Wb (^^ by a Gaussian with 
a width d, i.e., 

I^.(x.) - £^ /rf-^^/^'^^e(-')V2^\ (54) 
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With equation (^) this gives 

Plugging equation (153) into dSQ) we find with the same C as above that 



A 



3/2 



Now T„j has its maximum where 

In other words, given exactly the same functional form for both the longitudinal bunch 
form and the transverse bunch potential the delay Tm has its maximum at exactly 
A = d. This result is not surprising as now both ingredients of the longitudinal bunch 
potential enter in a symmetric fashion, i.e., with the same functional form and on equal 
footing. Then the symmetric case of their width parameters being equal is clearly a 
special point; here the delay is maximal. 

The slight offset to unity of the result obtained using the softcore potential can 
consequently be attributed to the fact that there two different functional forms were 
used for the longitudinal bunch profile and the transverse potential. Actually, these 
two quantities have different physically origins. 

This result that a spherical bunch with A = d leads to the biggest delay can be 
understood in the following way: when the bunch is prolate, i.e., A > d, the potential 
is flatter and slows down the ions less, whereas in a shorter, oblate bunch the distance 
2xi, over which the delay is accumulated, shrinks faster than the increasing effect of 
the steeper potential barrier. Of course, somewhere in between these two extremal 
cases there has to be a maximum, apparently when the bunch is spherical at the 
turning point. 

One might argue that for very short bunches the potential barrier should finally 
become so high that the ions come to a stop on top of the barrier. Then the delay 
would be high, too. But it is not only the length A of the bunch, which determines 
the potential's form, but also its width d: the curvature muj'^ (|42|) at the potential top 



is limited by the wider of the two functions contributing to Wb (see equation (40), be 
this the Gaussian ion distribution or the softcore Coulomb interaction. 

With the more realistic softened Coulomb potential we obtained an optimal bunch 
form of A slightly less than d, but for simplicity we will use A = d in the following. As 
this paper focusses on the general behavior of the bunch and not on highly accurate 
numbers this small difference of a few percent is clearly negligible. 



Using A = d, i.e., a fixed we see from equation (51) that the maximum 

delay for an otherwise fixed trap configuration scales with the adjustable parameters 
and A as 

A^ N 

^ra^^^^- (58) 

The result obtained above, that the most stable configuration is a bunch which is 
spherical at the turning point, is of central importance. As shown in part I Q 
synchronization is an interplay between the dispersion of the trap and the delay that 
the ions experience when the bunch turns around at the turning point. The trap 
dispersion is fixed by the mechanical and electrical setup of the trap. Only the delay. 
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which is the other central ingredient of synchronization, can adapt. It is striking that 
in its most stable configuration the bunch has a spherical form at the turning point. 
The bunch observed in the central part of the trap is then the "projection" of this 
spherical bunch out of the mirrors. 

To arrive at equation (|5l|), which is independent of the ions' initial conditions 
and is "local in the turning point" , it was important that both the top of the bunch 
potential and the ions' energy relative to the CM had the same quadratic dependence 
on the distance from the bunch's center. This is a consequence of the constant slopes 
of the mirror fields. But if the mirror is not linear, the bunch's potential at the turning 



point can still be expanded as in equation (41). If the bunch is not contracted linearly 



from Lij into A in the real experiment's mirror field, then the energy €{xi) of the ions 
will not be quadratic in Xi, but in any case it is a symmetric function of the distance 
from the CM. It can therefore at least be approximated by a quadratic dependence. 
Then Tm depends on Xi, but the delay is still maximal for A d. 

3.5. The bunch length in the central part of the trap 

With the above result, i.e., by setting yl = d, we can use the transformation ( pl| ) to 
predict the observed bunch length Lt as 



q \ dX 



x„ 



With the model potential of equation (|l|) we have (— f^) = F and, by expressing 
the ions' energy through the acceleration voltage U as Eq = qU, the observed bunch 
length becomes 



QTT 

— d = VSXrnrf . (60) 
r 

For the definition of Xm see ( pj] ) . In our model trap Lb is determined solely by the 
geometric trap parameters at the bunch's turning point: the penetration depth Xm of 
the ions into the mirror and the beam width d at this point. There a spherical bunch 
is formed, which is then "projected out" into the trap's central region. The position 
Xm of the turning point itself is of course determined by the ions' energy. 

Note that the observed bunch length does not depend on either the mass m of 
the ions or on their number A'', though the time delay Tm does (see equations ( pT] ) 
and (|5^)): the reason for this unexpected behavior is that the condition for maximal 
delay is that the integral in the derivative ( p3| ) be zero — a condition which only 
depends on the geometric properties d and A at the turning point. More or heavier 
ions lead to a higher delay, but do not change the bunch length of maximal delay. 

That the turning point is of central importance for the bunching could already 
be seen from the transformation ( pT| ) : the form of the potential between the field free 
region and the turning point has no effect on the relation between A and Lf,; what 
matters is the slope at the turning point. 

With the numerical values of the experiment, i.e., Eq = 4.2 keV, g = +1 and 
an average F = 80 kV/m we get Xm = 5.25 cm. In reference a beam radius of 
d ~ 0.5 mm was derived from trajectory calculations. These trajectory calculations 
also had predicted a spherical bunch at the turning point. With this value of d we 
predict an observed bunch length in the central part of the trap of Lb = 1.45 cm. Ar+ 
ions with the above Eq = 4.2 keV have a velocity of vq = 14.2-^, our Lb therefore 
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corresponds to a length in the time domain of Wi, = 0.1 fis. For comparison Pedcrsen 
et al measured a value of Wf, ~ 150 ns (figure 9(a) of reference ||] for Vi > 4.5 kV). 

The greatest uncertainty in comparing our results to the measurements stems 
from the greatly simplified form of our model potential. In the experiment only the 
voltage on the last of the five electrodes was varied, we can therefore not expect that 
the resulting field is linear. Due to geometrical limitations, i.e., the length L of the 
trap is fixed and the maximal Xm is limited in the experiment, it will be difficult to 



check the validity of equations (B9h or (60) over a wide range of parameters. Also the 



beam width d and the slope of the mirror potential {—§y) turning point have 

to be determined from trajectory simulations for each field configuration. 



With (y_l|) and (|13|) the kinematic regime a < 1 in which the trap has to be 



operated is described by the following condition (sec also equation (yi)| 

E,>!df or X,„^^>^ (61) 

The second form is a requirement on the trap's geometry: the mirrors have to be 
"long" enough, otherwise bunching cannot be observed. 

From this condition — which is another way of expressing the bunching criterion 
a < 1 — we see that there is a minimal bunch length determined by the length L of 
the field free region and the beam diameter 2c? at the turning point: 



Note that though the delay increases with the number of ions in the bunch (58), the 
observed itself is independent of N in this approximation and is determined by 
external parameters only. 



4. Numerical illustration of the behavior of the time delay 

In Section ^ we set up and determined with the quadratic approximation to 
the bunch potential. The result is described by equations ( ^ and (^. Then we 
linearized T ( ^5|) and found our central result that the time delay is maximal for 
A = d. Now we go back to equations ( p7|) and ([l5| ) and numerically calculate the 
time delay from them. This will both confirm our findings with the linearized T and 
illustrate the behavior of the bunch with respect to the most important parameters. 

The delay depends on a number of parameters, some of which are "externally 
given", i.e., fixed for the given experimental setup, while others describe how the 
bunch "reacts" to this environment. For instance, the size L of the trap, the energy 
Eq, the charge q and the mass m of the ions and the slope of the mirror fields F 
are fixed. Derived from those are the momentum pq = \/2mEQ, the time Tm ( pi)[ ) 
and the dispersion parameter a (p^), which we consider fixed, too. In the following 
these quantities will be set to their respective experimental values of Eq = 4.2 keV, 
F = 80 kV/m, q = +1 e and m = 40 amu, which were used in the initial experimental 
discovery of the bunching effect ||l|. With these values we get Tm = 1-48 ^s and 
a = 0.956. 

On the other hand we have the number of ions in the bunch A'', the observed 
length Lb and the length A in the mirror. There is no constraint on these quantities, 
so the bunch will try to adjust them to achieve the most stable configuration possible 
in the given trap regime. The two lengths, Lt and A, arc not independent, though. 
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Figure 3. Time delay caleulated from{|^), @ and (^) for a constant number 

of ions A' = 10^ and fixed beam diameters oi d = 0.2 mm ( ), 0.5 mm ( ) 

and 1 mm ( ). The arrows indicate for each curve where the bunch length is 

equal to the beam diameter, i.e., A = d. The thin broken curve plots for d = 0.5 
mm the weighted average of the delay from the equations of motion with the full 
bunch potential, i.e., without the harmonic approximation (see text). 

They are connected to each other by the transformation (^l]), determined by the slope 
of the mirror potential at the turning point. We view the minimal length in the mirror, 
A, as the more fundamental of these quantities, characterizing the minimal width in 
the mirror region where the ions cross the bunch; the experimentally observable Lh is 
merely the projection of this contracted length out of the mirror. We do not expect 
the absolute values which we predict for the observed Lb to be exact, but we expect 
that the important dependencies, governed by A, will be reproduced. 

The parameters we mentioned up to now are easily classified as fixed or variable. 
The beam radius d has a special role: it belongs to the fixed quantities, as it 
is determined by the electrostatic potentials on the mirror electrodes, but in our 
treatment it appears as a constituent of the bunch potential Wb (^o|) and it is closely 
intertwined with A, as we have seen above. We already know from equation (|5^) that 
the delay is maximal for A d. Nevertheless we will treat cZ as a variable parameter in 
the following numerical illustrations — just to see our results of the previous section 
confirmed. 

We will now in turn vary each of these three parameters. A, d and N , while the 
other two remain fixed. This will give three different views on the behavior of the 
bunch. 



4.I. Time delay vs. bunch length A and validity of the approximations 

First we look at the time delay as a function of the bunch length A for fixed N and 
d. We evaluate equations (|7|), ( fls] ) and (^) with iV = 10^ for three values of d — 
0.2 mm, 0.5 mm and 1 mm — and plot the results in figure ^ For short bunches, i.e., 
small A, and for long bunches the delay is small with a maximum in between. This 
maximum of r,„ occurs around A = d, where the length and the width of the bunch 
are equal. This is exactly the behavior predicted in the previous section. 
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Figure 4. Time delay Tm{xi) computed from tlie equations of motion (solid 
line) and the ion density p{xi) of the bunch at the turning point (broken curve) 
against the initial distance Xi (|33[), both for A = o! = 1 mm (cf. equation (^)). 
The solid horizontal line denotes the weighted average (r) of Tm, while the two 
vertical broken lines at Xi = ±A mark the width of the ion bunch. 



The linearized result (|51|), from which we derived the condition ^ = d, is 
indistinguishable from the full solution in figure]^; it differs by less than one percent for 
small A and coincides for larger A. This highlights the fact that with these parameters 
the delay of the ions through the bunch is a small perturbation compared to the 
momentum kick from the mirror potential. 

To check the validity of our harmonic approximation of the bunch potential (^Tj) 
we numerically evaluated Tm from equation (pj). This means that we still keep the 



modified three-stage model of section 3.1, where the interaction between the test ion 
and the bunch is confined to the second stage, during which the bunch is static. The 
position dependent momentum p{z) in equation (HJ) is calculated from its energy from 
the momentum kick €{xi) (^) and the bunch potential Wb{xi) of equation (^ ) as 



p{z) = ^y2E{z)m with E{z) = e{x,) - {Wb{0) - Wb{z)) . (63) 

E(z) is the energy of the ion above the bunch potential. The resulting delay 
Tm = Tm,{xi) uow dcpcnds on the ion's position in the bunch Xi (^) at the end 
of the first stage, during which the bunch was compressed. 

In figure ^ the delay r™ is plotted against the starting position Xi for yl = d = 1 
mm. For comparison also the contracted longitudinal bunch profile p{xi) is shown. We 
see that Tm decays much slower with increasing distance Xi from the bunch's center 
than p{xi): at Xi = 2.5 mm, where the bunch density becomes negligible, t„i is still 
half of the maximal value close to the center. Thus most of the ions have a t,„ of 
about the same magnitude. Consequently, the density weighted average 

(r) = j dxi p{xi) Tm{xi) (64) 

is a reasonable approximation over the whole length of the bunch. Also for different 
values of d the ratio between the widths of p and r remains unchanged, figure ^ is a 
representative example. 

The density averaged delay (t) for d = 0.5 mm is given in figure ^ as a thin broken 
curve. It is quite close to the result with the harmonic approximation to the bunch 
potential. Thus, both the approximation to Wb (fll[) and the subsequent linearization 



in the 


"Zajfman trap" II 




1 






= 1 1 1 1 1 IV 1 1 1 1 1 

- A = 0.5 mrrT 5v 


1 1 1 1 1 1 1 1 = 


0.1 


= A = 2 mm — 


- 

. A=d - 


U.U 1 


^A = 8mm ~ — ■ 




0.001 


1 


1 N 




0.1 

d [mm] 


1 10 



19 



Figure 5. Time delay as a function of tlie beam diameter for = 10^ and various 

fixed bunch lengtlis of yl = 0.5 mm ( ), 2 mm ( ) and 8 mm ( ). 

The enveloping solid line marks the time delay for A = d according to (|58[); the 
dots indicate the diameters, which correspond to the three bunch lengths shown. 
These correspond to values of Li, of 14, 29 and 58 mm, respectively. 



of T ((|48[) and (|49|)) are justified. For the following illustrations we will therefore 
continue to evaluate equations (37), (|45|) and (EgI). 



4-2. Time delay vs. beam diameter d 

When we now keep the bunch length fixed and vary the beam diameter the maximum 
of the delay for a spherical bunch shows up again. In figure ^ the delay is plotted 
for three different bunch lengths. It decreases with increasing beam diameter, but 
for fixed d the maximum is at A = d. This configuration is indicated by the solid 
enveloping line with a slope of —3/2 according to equation (^8|). 

From this plot it can be seen, too, that when the spatial dispersion requires a 
minimal delay of, e.g., ATi/A ~ 0.025 ns (see (^8|)) then the beam has to be focussed 
to a radius of no more than 1 mm. 

The necessary radial focussing can also be determined by inserting A ~ d and 



Arm > ATi (28) into equation (p l\) . The beam radius in millimeters has to be smaller 
than 



d[mm] < 5.02 x 10"^^^-^^ ^rj^ , (65) 

F AT^ '-^ 

when the mass m is given in atomic mass units (amu), the charge q in elementary 
charges, the slope F in kV/m and the dispersion ATi hi nanoseconds. 



4-3. Time delay vs. ion number N 

Another view of the behavior of the time delay is presented in figure |^. Now for given 
N and d the bunch length A is optimized to maximize the time delay. Panel (a) plots 
the maximal delay according to the optimized A. Panel (b) gives the corresponding 
observed bunch length if,. Both the maximal t,„ and Lf, are plotted vs. N for two 
different values of d. 

For small N the delay grows linearly with iV, as already derived in equation (p8|). 
This behavior is depicted by the thin solid lines of slope 1 in panel (a). In this linear 
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Figure 6. (a) Maximal time delay achievable at a given beam diameter and (b) 
the observed bunch length Li,, for which this maximal Tm is reached, as a function 
of the number of ions A'^. The thin solid lines indicate in (a) the linear relation 
between Tm and A'^jfcdland in (b) the constant A = d according to the linearized 
solution of section B.3.3|. 



regime the bunch is spherical, as can be seen in panel (b): Li, is constant, which 
implies that A is constant, too (see (p2[)). The thin solid lines in panel (b) mark the 
values of Lb, that result from inserting the two beam radii d = 0.5 mm and d ~ 0.2 
mm into (|32|). 

For higher ion numbers, however, the maximal achievable delay increases (panel 
(a)) and is reached for shorter and shorter bunches (panel (b)). Finally, the ion 
number, and therefore the density of the bunch, will become so high that the ions 
cannot pass over the bunch any more with the energy from the momentum kick. The 
ions then bounce back from the bunch and do not reach the other side any more. 
This is where 4po < Ltmuj, consequently our formulation for the crossing time T ( ^5|) 
breaks down. 

Inserting equations (|46|), ( |3^ ) and ( [l(i| ) into this condition 4po > Lbmuj for ( [l5| ) 
to be valid we see that this determines the maximal number iVmax as 

4V2^ d^F 

^max = ^TTTT ■ 66) 

Note that iVmax is proportional to the beam area via d^: this is consistent with our 
previous explanation that there is a maximal potential barrier from the contracted 
bunch that the ions can overcome with the energy from the momentum kick. 

Consequently there are — for a given d — two limits for the number of ions in a 
bunch: the lower limit is determined by the minimal delay necessary to compensate 
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for the spatial dispersion AT!; and the maximal stems from the requirement that 
the ions have to be able to pass over the contracted bunch's potential. 

Equation (|6^) in the previous section can be solved for the minimal number of 
ions necessary to support a given spatial dispersion AT^: 

Ntnin = -^r— ATi (67) 

1(1) ^rnq 

With the numbers from the experiment we calculate from ( |66|) and ( |67| ) that for a 
beam radius of c? = 0.5 mm and a dispersion of AT; =0.1 ns |l[| the number of ions 
may vary in the range of = 3 x 10^ . . .2 x 10^, i.e., over nearly three orders of 
magnitude. Within this range the bunch length is independent of the number of ions. 

This is consistent with the experimental observation that the bunch length is 
essentially independent of the density of the injected ion beam over at least three 
orders of magnitude Our estimate is more restrictive, especially A^max is bigger 
in the experiment. The reason lies in our harmonic approximation for Tm- as we see 
from figure |^ the delay is smaller for ions in the wings of the bunch, which means that 
for them the momentum kick increases faster than the height of the potential barrier. 
These ions will be able to pass over a bunch, which already reflects the ions close to 
its center. These will move away from the center of the bunch. Consequently the form 
of the bunch changes, allowing the bunch to be stable for higher ion densities at the 
turning point than we estimated here. In order to find the "true" A^max a dynamical 
model with a truly self consistent bunch potential has to be used; our fixed Gaussian 
form of the density is too simplified to give accurate numbers in this extreme regime. 

5. Momentum spread of the bunch 



In section 2.4 we stated that the bunch is heated up to a specific energy spread in the 
trap by the momentum kicks from the trap potential, no matter how small the initial 
energy spread at injection time had been. 

In "part I" we derived a relation between the maximal distance between the two 
ions a;niax and their maximal relative momentum Pmax in the central part of the trap. 



As explained in section 2.4 in the real trap the effective delay r^^ = t„i — ATi/ 4 has 
to be used, and therefore equation (37) of reference becomes 

p_, = ^ with /3=^./^^^^%:^-(l-«)2. (68) 

This relation is valid for a constant r^j. Our harmonic approximation to the bunch 
potential leads to the same constant delay for all ions of the bunch; consequently we 
can use this relation to relate the longitudinal momentum spread of the bunch Ap to 
its extension in the central part of the trap Lh as Ap = L^/p. 

With the parameters of the original experiment, i.e., Tm — 1.48 /is, a ~ 0.956 
and < 0.1 ns, the second term in the square root in the above equation ( |68| ) is 
smaller than the first by six orders of magnitude and can therefore be neglected. We 
see then, that Ap scales with the square root of the effective time delay t,'„ and is 
proportional to the bunch length L^: 



Ap^—J ^^ifc (69) 

^ T„VT„(l-a) ^ ^ 



Synchronization in the "Zajfman trap" II 



22 



This form highhghts that there are two different routes from the stable regime to the 
stability limit. One route is to increase the slope of the mirror potential to change 
the dispersion of the trap so that a — > 1. In this case the momentum spread grows 
only slowly as long as the trap is "far" in the stable regime, but as the mirror slope F 
approaches the maximal value Fmax = (equation ( p^ with a = I) Ap increases 
sharply and the bunch becomes drastically "hotter" . Right at the stability limit a ~ 1 
we get Ap ~ oo. The bunch is "blown apart" by its internal energy. Correspondingly 
the notion of a "momentum spread of the bunch" loses its meaning, because the 
bunch itself ceases to exist. This "heating" of the bunch is a collective effect of all 
ions together. 

The other route, which was discussed less extensively, is to relax the radial 
focussing of the beam to allow for stronger off-axis motion of the ions. This increases 
ATi, and therefore decreases t^j. Then Ap decreases with ATi/A T,n, i.e., the 
bunch becomes "cooler" closer to the stability limit. This happens because a stronger 
off-axis motion increases the spread of the periods of the individual ions through the 
trap, which effectively stretches the bunch in time and consequently decreases the 
interaction between the ions at the turning point. In contrast to the first collective 
route to instability via the longitudinal trap dispersion this second effect stems from 
the transverse dispersion due to different off-axis orbits of the individual ions. On 
the first route the bunch is "blown apart" because the bunching mechanism cannot 
confine the internal energy any more, while on the second route the bunch "boils off" 
ions that stray away too far from the bunch. Only those ions that have similar orbits 
through the trap stay together, thus reducing the internal momentum spread of the 
remaining bunch. 

We once again can insert the numbers from the original experiment, i.e., Tm ~ 
1.48 fis, F ~ 80 kV/m, q = +1 and a ~ 0.956 and find from figure || that 
Tm = 0.07 ns at d = 0.5 mm for N = 10^ ions. With a realistic ATi = 0.18 ns 
[01 and Xniax ^ Lb = 1-44 cm we calculate a momentum spread of Ap = 13.5 au. 
With m = 40 amu this translates into an energy spread in the laboratory frame of 
AE = J^Ap = 0.87 au = 24 eV at a beam energy of Ea = 4.2 keV. 



6. Comparing to the empirical criteria for bunching 

In reference three empirical criteria were given by Pedersen et al for synchronization 
to occur. We now show how these criteria are related to the results derived from our 
mapping approach about how and when synchronization occurs and about the size 
and behavior of the bunch. 

(i) The first criterion is the so called "kinematical criterion" which in ^ is given 
as ^ > 0. This criterion is obviously the same as our stability criterion of a < 1 for 
Tm > 0, see equations ( |l6| ) and (p5|). Both expressions state that the dispersion of the 
trap has to be positive, i.e., that a faster ion has a longer period. 

(ii) The second criterion of Pedersen et al , the " focussing criterion" , states that 
the interaction of the ions is important only in the mirror regions and that the beam 
diameter has to be compressed to about the same width as its length. We have shown 
in part I (t) that collisions in the central part are far more unlikely than collisions 
in the mirror and that only collisions in the mirror can synchronize the ions. We 
found the requirement of comparable longitudinal and transverse dimensions at the 
turning point when we estimated the bunch length in sections ^ and ^ a spherical 
bunch is most efficient in delaying the ions and therefore most stable against external 
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perturbations. With this rcsuh, i.e., A ^ d, and the transformation (3l]) between A 
and Lh we then derived an estimate for the observed bunch length (|9), which agrees 
with the prediction of Pedersen et al (see equation (29) of reference |S ]) . 

(iii) The third criterion of Pedersen et al , the "coUision criterion", states "[. . .] 
the colhsion probabihty at the turning point must ensure that the ions indeed lock their 
motion [...]" and "[...] too few collisions and too many lead to diffusion [...]". In 
our treatment the "number of collisions" corresponds to the magnitude of the effective 
time delay = r„i — > (see equation (pSh). As discussed in sections 4.2 and 4.3 



the delay t„i due to the collisions has to be bigger than the intrinsic dispersion ATi/4 
for synchronization to survive this perturbation, corresponding to the observation that 
too few collisions lead to diffusion. We want to emphasize again that the time spread 
AT„ of reference due to the ions' different energies is not a perturbation but a part 
of the ion dynamics. 

The upper limit of too many collisions, was illustrated in figure ^: if the bunch 
contains too many ions or is focussed too tightly, the ions bounce off the bunch 
potential, experience a negative delay and diffuse out of the bunch. In the framework 
of our approximations this upper limit is described by equation (Bq). 



7. Summary and conclusions 

In this paper we extended a recently proposed description of the synchronization 
effect in the "Zajfman trap" from describing the motion of two identical ions to 
the behavior of a macroscopic ion bunch: the motion of the two identical ions through 
the trap field was described in their CM system by three mappings, one for each 
of the different parts of the trap. When the paths of the ions cross in the mirror 
their repulsive interaction slows them down temporarily and the resulting time delay, 
together with the dispersion of the trap, is responsible for coupling their motion. Here 
we generalize this two ion approach to describe the motion of a (representative) test 
ion in a macroscopic bunch. The same stability condition applies as in the two ion 
case. The interaction between the ions is now specified through the size and form of 
the bunch. Together with a relation between the observed bunch length in the central 
part of the trap and its size at the turning point in the mirror the time delay can be 
calculated explicitly. 

We show that the time delay is maximal for a bunch which is spherical at the 
turning point and we calculate how the observed bunch length derived from this most 
stable configuration depends on the various parameters of the system. 

Not only docs our prediction of the observed bunch length match the experimental 
results within the level of approximation of our model, but we also confirm that the 
bunch length is essentially independent of the number of ions. The lower and the upper 
limits to this number are identified: the lower limit is connected to the spread of the 
ions' periods due to off-axis motion, as the time delay that the ions experience when 
they cross the bunch has to be big enough to compensate for it. At the upper limit 
the repulsive bunch potential of the compressed bunch at the turning point becomes 
too high and the ions cannot cross the bunch any more. Then they bounce back from 
the bunch, the resulting time delay changes its sign, the ions do not synchronize any 
more and consequently the bunch loses ions. 

The numerical evaluation of the time delay for various parameters illustrates the 
behavior of the bunch and confirms the validity of the approximations employed in 
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of our description. We make several new predictions that we hope wiU be tested 
experimentally: 

(i) To directly observe our central result of A = d at the turning point some kind 
of imaging technique has to be used. This would allow one to capture the form 
of the bunch at that moment when it comes to rest without having to resort to 
trajectory calculations. 

(ii) Much easier than the direct visualization of the turning bunch is to check the 
validity of equation (60), which is based on A ~ d. For a constant mirror slope 
F at the turning point we predict that the bunch length observed in the central 
part of the trap Lb scales with the square root of the beam diameter 2d at the 
turning point. For a fixed d, Lb would be inversely proportional to the square 
root of the slope F. However, Lb cannot become shorter than a minimal value, 
which depends on the total length L of the field free region of the trap according 
to (H). 

(lit) We found that the bunch length is determined by the mirror slope at the turning 
point only and not by the dispersion, see (|59|). This can be verified with a 
nonlinear mirror potential, which keeps the slope at the turning point constant 
for different trap dispersions. 

(iv) The observed bunch length is independent of the number of ions in the bunch 
within certain limits, which are given by equations ( |66| ) and (^). The two limits 
scale differently with the parameters of the trap and the ions. 

(v) With equation ( |69| ) we derived a relation between the bunch size and its 
longitudinal momentum spread. According to this equation there are two ways 
that a bunch can become unstable, either via the longitudinal or via the transverse 
dispersion of the trap. The two routes, depending on the mirror slope and the 
beam diameter, respectively, have opposite effects on the momentum spread. 

Finally we explain the empirical conditions for synchronization given by Pedersen 
et al in the framework of our model. 

In this paper we considered only a static bunch, i.e., the bunch exists initially 
and is either stable or not. Future work therefore has to focus on the dynamical 
properties of the synchronization effect, as, e.g., how and at which rate the bunch 
acquires additional ions from which parts of phase space. Other open questions are 
related to the issues of ions of different masses: what is the mass resolution of a 
"Zajfman trap" operated as a mass spectrometer? Do bunches of different ions "stick 
together" or can one have more than one bunch in the trap independently? 

Another set of questions is connected to the internal dynamics of the bunch 
and whether there is an application for this "inner trap" moving through the outer 
macroscopic trap. 

It is not clear at the moment if our simple model can be used to answer these 
questions or if a more detailed description of the many-ion dynamics has to be used. 
Nevertheless, our simple model can explain the underlying mechanism and the stability 
criteria and provide insight into how the macroscopic bunch reacts to the externally 
given conditions and parameters. 
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